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. Recently, a model f(R) theory is proposed [1] which is cosmologically viable and 

distinguishable from ACDM. We use chameleon mechanism to investigate viability of the 
q ■ model in terms of Solar System experiments. 
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O 1 Introduction 

On ' 

There are strong observational evidences that the expansion of the universe is accelerating. 
These observations are based on type la supernova [2] , cosmic microwave background radiation 
[3], large scale structure formation [4], weak lensing [5], etc. The standard explanation invokes 
an unknown component, usually referred to as dark energy. It contributes to energy density 
of the universe with Q d = 0.7 where fid is the corresponding density parameter, see e.g., [6] 
and references therein. The simplest dark energy scenario which seems to be both natural 
and consistent with observations is the ACDM model in which dark energy is identified as a 
cosmological constant [6] [7] [8]. However, in order to avoid theoretical problems [7], other 
scenarios have been investigated. Among different scenarios there are modified gravity models 
[9], in which one modifies the laws of gravity whereby a late time acceleration is produced 
without recourse to a dark energy component. One family of these modified gravity models is 
obtained by replacing the Ricci scalar R in the usual Einstein-Hilbert Lagrangian density for 
some function f(R). These models are cosmologically acceptable if they meet some certain 
conditions. The most important ones is that they should follow a usual matter-dominated era 
preceding a late-time accelerated stage. Several models are proposed that admit cosmological 
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solutions with accelerated expansion of the universe at late times [10]. However, among all 
cosmologically viable f(R) theories there is still an important issue to be pursued, they must 
be probed at Solar System scale. In fact, changing gravity Lagrangian has consequences not 
only in cosmological scales but also in galactic and Solar System ones so that it seems to be 
necessary to investigate the low energy limit of such f(R) theories. 

Early works on the weak field limit of f(R) theories led to negative results. Using the equiv- 
alence of f(R) and scalar-tensor theories [11] [12] [13], it is originally suggested that all f(R) 
theories should be ruled out [14] since they violate the weak field constraints coming from 
Solar System experiments. This claim was based on the fact that f(R) theories (in the metric 
formalism) are equivalent to Brans-Dicke theory with uo = while observations set the con- 
straint uj > 40000 [15]. In this case the post-Newtonian parameter satisfies 7 = \ instead of 
being equal to unity as required by observations. Later, it was noted by many authors that 
for scalar fields with sufficiently large mass it is possible to drive 7 close to unity even for null 
Brans-Dicke parameter. In this case the scalar field becomes short-ranged and has no effect 
at Solar System scales. Recently, it is shown that there exists an important possibility that 
the effective mass of the scalar field be scale dependent [16]. In this chameleon mechanism, 
the scalar field may acquire a large effective mass in Solar System scale so that it hides local 
experiments while at cosmological scales it is effectively light and may provide an appropriate 
cosmological behavior. 

In the present work we intend to use this criterion to investigate constraints set by local ex- 
periments on a cosmologically viable model recently proposed by Miranda et al [1]. It is a two 
parameter f(R) model which is introduced in the form 

f(R) = R-aR 1 \n(l + ^) (1) 

where a and R\ are positive parameters. This model is reduced to general relativity for a = 0. 
A cosmologically viable f(R) model must start with a radiation-dominated era and have a 
saddle point matter-dominated phase followed by an accelerated epoch as a final attractor. 
This is formally stated by introducing the parameters m = R d f J£ R and r = -R^j^ [17]. 
The cosmological dynamics of f(R) models can then be understood by considering m(r) curves 
in the (r, m) plane. Using this criteria, the authors of [1] showed that the model (1) satisfies all 
the cosmological requirements for a > 1 and regardless of R\ t Here we will focus on viability 
of the model in terms of local gravity experiments. We will show that these local experiments 
rule out the model as an explanation for the current accelerated expansion of the universe. 



tin the model proposed in [1], it is stated that the relation (1) is a special case of a general parametrization 
unifying the models of [18] and [19]. In the latter models there is a parameter with the same dimension of R\ 
which is taken to be of the same order of the presently observed cosmological constant. The relation between 
this parameter and R\ and also the relevance of R\ with the time of the beginning of the acceleration phase 
are important issues that are not explicitly addressed in [1]. 
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2 Chameleon Mechanism 



In this section we offer a brief review of the chameleon mechanism. We consider the following 
action* 

S = I J d 4 xV=9 f(R) + S m (g^, ij) (2) 

where g is the determinant of g^ u , f(R) is an unknown function of the scalar curvature R and 
S m is the matter action depending on the metric g^ v and some matter field ip. We may use a 
new set of variables 

9\xv = V 9nu (3) 

0=-^lnp (4) 

where p = ^ = f'(R) and (3 = ^/|. This is indeed a conformal transformation which 
transforms the above action in the Jordan frame to the Einstein frame [11] [12] [13] 

X - J d 4 xy/=g- {R - g^QA dA ~ 2^(0)} + S m (9^e 2/3<t> , ij>) (5) 
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In the Einstein frame, is a scalar field with a self-interacting potential which is given by 

V{4>) = \e- 2P<t> {M<t>)\ ~ e~ W /(rb(0)])} (6) 

where r(p) is a solution of the equation /' [r(p)]—p = [11]. Note that conformal transformation 
induces the coupling of the scalar field with the matter sector. The strength of this coupling 
f3, is fixed to be and is the same for all types of matter fields. In the case of such a 
strong matter coupling, the role of the potential of the scalar field is important for consistency 
with local gravity experiments. When the potential satisfies certain conditions it is possible 
to attribute an effective mass to the scalar field which has a strong dependence on ambient 
density of matter. A theory in which such a dependence is realized is said to be a chameleon 
theory [16]. In such a theory the scalar field can be heavy enough in the environment of 
the laboratory tests so that the local gravity constraints suppressed even if f3 is of the order 
of unity. Meanwhile, it can be light enough in the low-density cosmological environment to be 
considered as a candidate for dark energy. 

Variation of the action (2) with respect to g^ v and 0, gives the field equations 

GV = <9 M d u <p - ^g^d^ <9 7 - V((p)g^ + % v (7) 

dV 

□0 - £l = -FT (8) 

where 

f = (g) 



*We use the unit (%-kG)- 1 = 1. 
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and T = g^T^. Covariant differentiation of (7) and the Bianchi identities give 

V%, = (3f d u <p (10) 

which implies that the matter field is not generally conserved and feels a new force due to 
gradient of the scalar field. Let us consider T^ v as the stress-tensor of dust with energy density 
p in the Einstein frame. In a static and spherically symmetric spacetime the equation (8) gives 

d 2 <P | 2d<f> _ dV eff (<f>) 
df 2 f df d(j) 

where f is distance from center of the symmetry in the Einstein frame and 

Veff(<t>) = V(4) - \pe-^ (12) 

Here we have used the relation p = e^^p that relates the energy densities in the Jordan 
and the Einstein frames. We consider a spherically symmetric body with a radius f c and a 
constant energy density p in (f < f c ). We also assume that the energy density outside the body 
(f > f c ) is given by p out . We will denote by ip in and ip out the field values at two minima of 
the effective potential V e //(0) inside and outside the object, respectively. They must clearly 
satisfy V e ^{ifi n ) = and V e ^{ip ou t) = where prime indicates differentiation of V e //(0) with 
respect to the argument. As usual, masses of small fluctuations about these minima are given 
by m in = [V e ff(ipi n )]2 and m out = [V e ff((p out )] 5 which depend on ambient matter density. A 
region with large mass density corresponds to a heavy mass field while regions with low mass 
density corresponds to a field with lighter mass. In this way it is possible for the mass field 
to take sufficiently large values near massive objects in the Solar System scale and to hide the 
local tests. For a spherically symmetric body there is a thin-shell condition 



Ar e _ Lp out — i^ in 

where $ c = M c /87rf c is the Newtonian potential at f = f c with M c being the mass of the body. 
In this case, equation (11) with some appropriate boundary conditions gives the field profile 
outside the object [16] 

(3 3Af c M c e- m °^-rc) 
0(r) = ~4^~ f + Vout ( } 



3 The Model 

The function f(R) in the Jordan frame is closely related to the potential function of the scalar 
degree of freedom of the theory in the Einstein frame. Any functional form for the potential 
function corresponds to a particular class of f(R) theories. To find a viable function f(R) 
passing Solar System tests one can equivalently work with its corresponding potential function 
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in the Einstein frame and put constraints on the relevant parameters via chameleon mechanism. 
Taking this as our criterion, we write potential function of the model (1) 



vw = 2^ e " w {« ln (r^w) - eW4> - ( a - !)> ( 15 ) 

Assuming that << 1, one can find the solution of V e 'jj(0) = by substituting (15) into (12) 
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^ = —{-C ± VC^AAB} (16) 
where 

A = 2f3 2 R 1 (l-2a) (17) 

B = —aR-y (18) 

C = P[p + R 1 (a-l)-2aR 1 lna] (19) 

In the following we shall consider thin-shell condition together with constraints set by equiva- 
lence principle and fifth force experiments. 

1. Thin — shell condition 

In the chameleon mechanism, the chameleon field is trapped inside large and massive bodies 
and its influence on the other bodies is only due to a thin-shell near the surface of the body. 
The criterion for this thin-shell condition is given by (13). If we combine (13) and (16) we 
obtain 

Arv 1 



fr 12$ C A 



{(pin - pout) ± V (pout + a) 2 ~ b =F V (pin + a) 2 - b } (20) 



where pi n and p ou t are energy densities inside and outside of the body in the Jordan frame. 
Here a = 4 — p and b = AaRl(2a — 1). In weak field approximation, the spherically symmetric 
metric in the Jordan frame is given by 

ds 2 = -[1 - 2X(r)]dt 2 + [1 + 2Y(r)]dr 2 + r 2 dQ 2 (21) 

where X(r) and Y(r) are some functions of r. There is a relation between r and f so that 
f = p l / 2 r. If we consider this relation under the assumption m out r < 1, namely that the 
Compton wavelength m~„\ is much larger than Solar System scales, then we have f ps r. In 
this case, the chameleon mechanism gives for the post-Newtonian parameter 7 [24] 

3 _ Aik 9 Ar 

7 = ^7-1-'— (22) 



3 + ^ 3 r c 

Tc 



We can now apply (20) on the Earth and obtain the condition that the Earth has a thin-shell. 
To do this, we assume that the Earth is a solid sphere of radius R e = 6.4 x 10 s cm and mean 
density p e ~ 10 gr/cm 3 . We also assume that the Earth is surrounded by an atmosphere with 
homogenous density p a ~ 10~ 3 gr/cm 3 and thickness 100km. For simplifying equation (20), 
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we proceed under the assumption p in , p out « R±a§. We will return to this issue later. In this 
case, equation (20) simplifies to 

A-R e ^ 1 (Pin - Pout) / 2 ox 

R e ~ 4$ e i? 1 (l-2a) { ' 

where $ e = 6.95 x 10~ 10 [25] is Newtonian potential on surface of the Earth. The tightest 
Solar System constraint on 7 comes from Cassini tracking which gives | 7 — 1 |< 2.3 x 10~ 5 
[15]. This together with (22) and (23) yields 

(l-2«) > 10 13 (^) (24) 
Ri 

With p in = p e = 7 x 10~ 28 cm" 2 , this is equivalent to i?i(l - 2a) > 1(T 15 cm~ 2 . 
2. Equivalence principle 

We now consider constraints coming from possible violation of weak equivalence principle. 
We assume that the Earth, together with its surrounding atmosphere, is an isolated body 
and neglect the effect of the other compact objects such as the Sun, the Moon and the other 
planets. Far away the Earth, matter density is modeled by a homogeneous gas with energy 
density pa ~ 10~ 24 gr / cm 3 . To proceed further, we first consider the condition that the atmo- 
sphere of the Earth satisfies the thin-shell condition [16]. If the atmosphere has a thin-shell 
the thickness of the shell (AR a ) must be clearly smaller than that of the atmosphere itself, 
namely AR a < R a , where R a is the outer radius of the atmosphere. If we take thickness of 
the shell equal to that of the atmosphere itself AR a ~ 10 2 km we obtain < 1.5 x 10~ 2 . 
It is then possible to relate = and = V Q^f a where tp e , tp a and ipo are the field 

values at the local minimum of the effective potential in the regions r < R e , R a > r > R e and 
r > R a respectively. Using the fact that newtonian potential inside a spherically symmetric 
object with mass density p is $ oc pR 2 , one can write $ e = 10 4 <3> a where <3> e and $ a are 
Newtonian potentials on the surface of the Earth and the atmosphere, respectively. This gives 
AR e /R e PS 10 -4 AR a /R a . With these results, the condition for the atmosphere to have a 
thin-shell is 

ARe < 1.5 x 1(T 6 (25) 



R 



e 



The tests of equivalence principle measure the difference of free-fall acceleration of the Moon 
and the Earth towards the Sun. The constraint on the difference of the two acceleration is 
given by [15] 



< 10- u (26) 

where a m and a e are acceleration of the Moon and the Earth respectively and is the 
Newtonian acceleration. The Sun and the Moon are all subject to the thin-shell condition [16] 
and the field profile outside the spheres are given by (14) with replacement of corresponding 



3 It can be easily checked that our main results do not change when p in , p out » Rxa. 
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quantities. The accelerations a m and a e are then given by [16] 

a e ^{l + 18/? 2 (^) 2 ^} (27) 
A R cE> 2 

a m » a„{l + 18/? 2 (^) 2 ^} (28) 

where $ e = 6.95 x lO" 10 , $ m = 3.14 x 10~ n and $ s = 2.12 x 1(T 6 are Newtonian potentials 
on the surfaces of the Earth, the Moon and the Sun, respectively [25] . This gives a difference 
of free-fall acceleration 



Combining this with (26) results in 

AR, 



= (0.13) (3 2 (=P) J (29) 



< 6.74 x 10~ b (30) 



which is of the same order of the condition (25) that the atmosphere has a thin-shell. Taking 
this as the constraint coming from violation of equivalence principle and combining with (23), 
we obtain 



(l-2a)>(10 14 )(^) (31) 
which is not much different from (24). 

3. Fifth force 

The potential energy associated with a fifth force interaction is parameterized by a Yukawa 
potential 

U(r) = 32 

where mi and m 2 are masses of the two test bodies separating by distance r, e is strength 
of the interaction and A is the range. Thus fifth force experiment constrains regions of (e, A) 
parameter space. These experiments are usually carried out in a vacuum chamber in which 
the range of the interaction inside it is of the order of the size of the chamber [16], namely 
A ~ R V ac- The tightest bound on the strength of the interaction is e < 10~ 3 [26]. Inside the 
chamber we consider two identical bodies with uniform densities p c , radii r c and masses m c . If 
the two bodies satisfy the thin-shell condition, their field profile outside the bodies are given 
by 

,, . (3 3 Ar c m c e ~ r/Rvac , 00 , 

W = — A 1" Vvac (33) 

Then the corresponding potential energy of the interaction is 

v(r) = _ 2/j2( 3^ )2 

The bound on the strength of the interaction translates into 

2/3 2 (^) 2 < 10- 3 (35) 
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One can write for each of the test bodies 



Arc _ 1 (Pc-Pvac) m) 
r c ~ 4$ CJ R!(l-2a) 1 ; 

where p vac is energy density of the vacuum inside the chamber. In the experiment carried 
out in [26], one used a typical test body with mass m c 40gr and radius r c ps 1cm. These 
correspond to p c ~ 9.5gr/cm 3 and <3> c ~ 10~ 27 . Moreover, the pressure in the vacuum chamber 
was reported to be 3 x 1CT 8 Torr which is equivalent to p vac ~ 4.8 x 10~ 14 gr/cm 3 . Substituting 
these into (36) and combining the result with (35) gives the bound 

(l-2a)>(10 27 )(^) (37) 
which is equivalent to R±(l — 2a) > 1CT 1 cm~ 2 . 



4 Discussion 

We have discussed viability of the f(R) model proposed in [1] in terms of local gravity con- 
straints. We have used the correspondence between a general f(R) theory with scalar field 
theories. In general, in the scalar field representation of a f(R) theory there is a strong cou- 
pling of the scalar field with the matter sector. We have considered the conditions that this 
coupling is suppressed by chameleon mechanism. We have found that in order that the model 
(1) be consistent both with fifth force and equivalence principle experiments, the two parame- 
ters a and R\ together should satisfy the condition Ri(l — 2a) > 1CT 1 cmr 2 . To have a bound 
on the parameter a, one should attribute a physical meaning to the dimensional quantity R\. 
Following the models proposed by Hu and Sawicki [18] and Starobinski [19], if we take it as the 
same order of the observed cosmological constant A obs ~ 1CT 58 cm -2 we obtain \2a — 1| > 10 57 . 
It should be pointed out that this result is obtained under the assumption p in , p out « R\a. To 
understand the relevance of this assumption, let us consider the case that pi n , p ou t ~ R±a. Tak- 
ing energy density inside the Earth as a typical energy density in the Solar System, we obtain 
R\a ~ 10~ 28 cm~ 2 . In the model (1), if the coefficient R\a is so small then it would be hardly 
distinguishable from general relativity. As the point of view of local gravity experiments, a 
viable f(R) model should simultaneously satisfy Solar System bounds as well as exhibit an 
appropriate deviation from general relativity. This requires that R±a » 1CT 28 cm -2 , which 
confirms the assumption that our results are based on. 

The last point we wish to remark is that, as reported by the authors of [1], the condition that 
the universe pass through a matter-dominated epoch and finally reach a late-time accelerated 
phase is that a > 1 regardless of R\. This seems not to be consistent with our results in the 
context of the chameleon mechanism. Although our analysis do not place any experimental 
bound on the parameter Ri, however for R\ > the relation (37) implies that a G (— oo |] 
which is out of the range reported in [1]. 
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